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The  paper  presents  an  initial  move  to  develop  a  scalable  high-degree-of-freedom  model  of  the  parafoil- 
payload  system.  The  intention  is  to  develop  the  tool  capable  of:  i)  determining  basic  system’s  geometry 
parameters  by  observing  the  video  data  of  the  real  descend,  ii)  readjusting  the  nominal  aerodynamic  and 
control  coefficients  incorporated  into  the  well-established  equations  of  motions,  and  iii)  performing  model 
identification  to  tune  numerous  relative  variables  to  achieve  the  best  fit  with  the  real  drop  data  if  available. 
Since  in  the  certain  way  such  a  tool  would  represent  some  kind  of  generalization  of  the  modeling  efforts 
undertaken  so  far,  the  present  paper  starts  from  a  comprehensive  review  of  publications  devoted  to  the 
modeling  of  parafoil-payload  systems.  The  paper  then  briefly  addressed  the  current  stage  of  the  development 
of  a  scalable  model.  In  anticipation  of  real  drop  data  to  validate  the  approach  paper  ends  with  conclusions. 


I.  Introduction 

DUE  to  their  ability  to  cover  large  horizontal  distances  from  the  drop  point  and  their  controllability,  gliding 
parachutes  offer  considerable  scope  for  military  application  in  the  delivery  of  troops,  stores  and  munitions, 
enhanced  safety  for  aircrew,  etc.  One  concept  in  the  gliding-parachute  category  which  has  received  attention  in  mid 
60’s  is  a  ram-air-inflated  fabric  wing  (called  the  Jalbert  load- distributing  flares  parafoil).  It  has  a  rectangular 
planform  and  an  airfoil  cross-section  with  an  opening  at  the  leading  edge  to  allow  ram  air  to  enter  and  inflate  the 
wing  to  the  desired  shape.  The  payload  is  suspended  on  lines  below  the  wing.  Some  of  the  earliest  projects  involving 
the  usage  of  parafoil-based  systems  were  discussed  by  Nikolaides  and  Knapp.1'7  Lingard8  provided  comprehensive 
bibliography  of  ram-air  literature  up  to  1980.  Some  definitions  that  in  a  quarter  of  a  century  have  become  sort  of 
standard,  discussion  on  modem  parafoil  parameters  and  their  design  were  presented  by  Poynter.9  For  the  past  40 
years  a  lot  of  models  of  different  parafoil-payload  systems  were  developed  to  address  the  variety  of  issues.  Lately, 
the  high-fidelity  8  and  9  degrees-of  freedom  (DoF)  models  were  also  created. 

This  paper  describes  current  status  of  an  attempt  to  develop  a  scalable  model  of  a  generic  parafoil-payload 
system  based  on  the  model  that  has  already  been  developed  earlier  for  the  Pegasus  payload  delivery  system  by  FXC 
Corp.10  In  what  follows  the  extensive  review  of  the  aerodynamic  data  available  for  different  parafoils  as  well  as 
numerous  attempts  to  model  parafoil-payload  systems  are  given  first.  Then  a  brief  description  of  the  software 
developed  so  far  is  presented.  This  tool  also  employs  the  system  identification  tools  to  better  match  the  real  drop 
data  if  available.  These  tools  however  are  considered  in  a  separate  paper.11 


II.  Historic  Background 


A.  Wind-Tunnel  Tests 

Since  mid  60 ’s  when  ram-air  parafoils  were  first  introduced  many  studies  have  been  carried  out  to  establish 
their  basic  aerodynamic  performance,  range  of  stability  and  optimal  design,  to  develop  relevant  models  and  recently 
control  laws  and  instmmentation  for  autonomous  payload  delivery  were  addressed.  The  design  of  ram-air  gliding 
parachute  canopies  employs  more  airplane  principles  than  parachute  technology  therefore  the  designers  employ  texts 
on  basic  aerodynamics.  One  problem  that  will  be  encountered  is  that  most  airplanes  do  not  fly  as  slowly  as  the  ram- 
air  canopy  (20-25 mph).  Consequently,  very  little  has  been  done  in  this  area  and  the  information  one  seeks  is  usually 
off  the  charts.  Therefore  in  what  follows  we  present  an  overview  of  different  efforts  on  modeling  (in  total,  over  a 
hundred  conference  papers,  technical  reports  and  journal  publications  have  been  collected  and  reviewed)  and  we 
start  from  the  review  of  experimental  data  available. 
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The  first  preliminary  wind-tunnel  tests  and  free-flight  tests  on  this  type  of  wing  were  carried  out  in  wind  tunnel 
of  university  of  Notre  Dame  by  Nikolaides  et  al.  in  1964. 1,2  Burk  et  al.  presented  an  investigation  of  the  low-speed 
static  aerodynamic  characteristics  of  three  ram-air- inflated  low-aspect-ratio  parafoils  (0.64,  0.78  and  0.94 
respectively)  in  the  NASA  Langley  Research  Center  full-scale  wind  tunnel.12  The  aerodynamic  coefficients  included 
lift  CL ,  drag  CD  and  side-force  CY  coefficients;  rolling  C, ,  pitching  Cm  and  yawing  Cn  moment  coefficients; 
effective-dihedral  Clf3 ,  directional  stability  Cn/3  and  side-force  parameters  CYJ3 .  The  moment  reference  center  was 

located  4.0  percent  of  the  wing  chord  aft  of  the  leading  edge  of  the  wing  and  123  percent  of  the  wing  chord  below 
the  lower  surface  of  the  wing.  The  results  indicated  that  the  untrimmed  maximum  lift-drag  ratios  of  the  wings 
ranged  from  1.9  to  2.7  and  that  the  untrimmed  maximum  lift  coefficients  varied  from  about  0.9  to  1.1.  The  results  of 
tests  of  two  of  the  wings  selected  for  a  study  of  their  stability  characteristics  indicated  that  they  were  statically 
longitudinally  stable  over  the  entire  test  angle-of-attack  range  of  0°  to  70°.  Both  of  these  wings  also  had  generally 
satisfactory  static  lateral  stability  characteristics  over  the  entire  test  angle-of-attack  range. 

Ware  et  al.  extended  previous  results  to  11  new  models  with  higher  aspect  ratios  from  1.0  to  3.0. 13  Lift-drag 
ratio  for  the  constant-area  series  was  reported  to  be  about  2.5  for  all  aspect  ratios  (about  the  same  as  for  a  full-scale 
version),  and  the  maximum  lift-drag  ratio  for  the  constant  chord  series  varied  from  2.2  to  3.3  with  increasing  aspect 
ratio.  It  was  also  noted  that  for  the  wing  alone  (without  suspension  lines  the  lift-drag  ration  was  in  average  of  one 
unit  higher).  The  data  which  contained  the  same  list  of  parameters  plus  profile-drag  Cm  and  induced-drag  CDi 
coefficients,  was  referred  to  the  wind  system  of  axes  and  based  on  the  laid-out-flat  wing  area  and  chord  of  span.  The 
moment  reference  was  the  confluence  point  of  the  suspension  lines,  and  the  one  used  in  stability  and  control  tests 
was  either  the  25 -percent  chord  of  the  lower  surface  of  the  wing  or  an  analytically  determined  center  of  gravity 
located  1.5  spans  below  the  wing  and  positioned  longitudinally  so  as  to  trim  the  configuration  at  the  maximum  value 
of  lift-drag  ratio. 

Numerous  parafoil  configurations  which  included  variations  in  aspect  ratio  (0.5  to  3.0),  airfoil  thickness  (0.126 
to  0.201),  airfoil  shape,  pennant  design,  leading  edge  opening,  control  surface  (flaps)  deflections  et  al.  were  tested 
by  Nikolaides.1,5  Also,  tests  were  carried  out  for  a  wide  range  of  parafoil  size,  wind  tunnel  velocity,  and  rigging.  The 
list  of  parameters  was  broadened  to  include  lift-curve  slope  CLa  ,  pitching  stability  parameter  Cma  ;  pitch  damping 
Cmq  and  lag  moment  Cd  coefficients;  drag  coefficient  of  lines  CDn  (which  cause  as  much  as  35%  of  the  parasitic 

drag).  Data  from  various  wind  tunnel  testing  programs  has  been  reduced  to  a  common  basis  and  numerous 
comparison  plots  illustrating  effects  of  aspect  ratio,  velocity,  trim  angle,  control  deflection  and  configuration  have 
been  presented  (see  also  Ref.2  where  the  results  of  the  real  drops  compared  to  that  of  experimental  prediction  are 
given). 

This  experimental  research  confirmed  that  parafoil  is  similar  to  airplane  wing.  It  has  positive  lift  at  zero  angle  of 
attack  and  shows  lift  down  to  about  -5°.  The  lift  curve  is  quite  linear  with  angle  of  attack.  Increasing  the  aspect  ratio 
increases  the  lift  curve  slope  and  improves  the  lift  to  drag  ratio.  The  parafoil  happens  to  have  static  and  dynamic 
stability  in  all  modes  of  flight  in  pitch,  yaw,  and  roll.  Because  of  its  configuration  and  flexibility,  it  does  not  exhibit 
the  stall  characteristics  of  the  aeroplane  wing  at  large  angles  of  attack.  Instead  the  lift  falls  off  gently  and  thus  the 
parafoil  may  also  be  safely  flown  at  very  large  angles  of  attack  (70°+). 

These  wind  tunnel  investigations  conducted  in  the  1960’s  in  the  University  of  Notre  Dame  2’x2’  test  section 
and  in  the  NASA  Langley  Research  Center  30’x60’  (elliptic)  test  section  were  conducted  on  models  at  relatively  low 
wing  loadings  (l-2lb/ft2)  and  small  size  models  up  to  300ft2.  They  constitute  the  majority  of  the  experimental  data 
available  for  small  (personnel  size)  systems.  Going  beyond  the  personnel  sized  canopies  (175  to  340ft2),  some  very 
limited  research  has  been  done  on  ram  air  inflation  systems  up  to  3,200 ft2. 

Since  the  ram-air  inflated  wing  has  many  proved  advantages  over  the  more  conventional  Parachute  system  such 
as  low  weight,  high  maneuverability  and  capability  to  flare  for  a  soft,  stable  landing,  in  late  80 ’s  Pioneer  Aerospace 
Corporation  chose  it  for  the  Advanced  Recovery  Systems  (ARS)  for  the  Next  Generation  Space  Transportation 
System.14,15  The  canopy  size  required  for  this  test  program  must  go  far  beyond  any  that  have  been  previously 
studied.  The  full  scale  prototype  (10,800 ft2,  5.56lb/ft2)  exceeds  the  size  of  3,200 ft2  by  338%.  Due  to  the  lack  of  data 
for  large-size  parafoils  a  large  scale  wind  tunnel  tests  in  the  NASA  Ames  Research  Center  80 ’x  120’  test  section  of 
the  National  Full-scale  Aerodynamic  Complex  at  Moffett  Field,  California  were  conducted  to  establish  a  data  base 
of  large  (1,200 ft2)  ram  air  inflated  wings.16  Two  models  were  tested  during  this  investigation  -  the  primary  test 
article,  a  1/9  geometric  scale  model  with  wing  area  of  1,200 ft2  and  secondary  test  article,  a  1/36  geometric  scale 
model  with  wing  area  of  300 ft2,  both  of  which  had  an  aspect  ratio  of  3.0.  The  test  results  showed  that  large  size 
models  were  statically  stable  about  a  model  reference  point  at  angles  of  attack  from  2  to  10  degrees.  The  maximum 
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lift-drag  ratio  varied  between  2.9  and  2.4  for  increasing  wing  loading.  (The  flight  test  results  of  the  developed  full- 
scale  X-38  system  were  discussed  in  details  nine  years  later  by  Iacomimi  et  al.,  Martin,  Madsen  and  Cerimele.17'23) 

B.  First  Modeling  Efforts 

To  this  point  all  aforementioned  research  provided  the  only  experimental  data  available  to  develop  different 
models.  Based  on  this  data  Lingard  applied  low  aspect  ratio  wing  theory  to  establish  a  model  from  which  the  glide 
performance  of  the  most  successful  gliding  parachute,  the  ram-air,  may  be  understood.24,8  His  predictions  were 
shown  to  compare  well  with  experiment.  Optimum  performance  for  a  conventional  ram-air  parachute  was  shown  to 
be  close  to  3:1,  the  value  obtained  in  practice.  Improvements  in  aerodynamic  efficiency  with  increasing  aspect  ratio 
were  shown  to  be  largely  negated  by  increased  line  drag.  Optimum  aspect  ratio  and  line  length  were  defined. 
Gliding  efficiency  was  shown  to  be  strongly  influenced  by  the  position  of  the  suspended  load  under  the  canopy  and 
the  ideal  location  was  discussed. 

Goodrick  used  a  3 -DoF  model  to  analyze  both  static  and  dynamic  longitudinal  stability  of  high-performance 
gliding  airdrop  systems  (parafoils).25  He  also  presented  results  of  a  6-DoF  simulation  and  its  comparison  with  the 
real  drop  (although  equations  of  motion  were  omitted).26,27 

Lingard  first  derived  equations  of  motion  (EoM)  of  gliding  parachute  system  by  resolving  the  forces 
horizontally  and  vertically  giving  relationships  for  total  horizontal  and  vertical  velocities.8  These  relationships  were 
to  be  numerically  integrated  to  calculate  two-dimensional  flight  trajectories.  It  accommodates  atmospheric  wind 
effect  also.  Lingard  used  this  model  to  study  the  effect  of  wind  on  range  covered  by  the  system.  Next,  two- 
dimensional  3 -DoF  dynamic  model  was  developed  which  was  defined  by  translational  velocity  components  and  the 
angular  pitching  component.  This  model  was  capable  of  simulating  effect  on  flight  dynamics  of  changes  in  canopy 
size,  line  length,  rigging  and  store  mass,  and  drag  along  with  changes  in  wind  conditions.  It  was  shown  that 
improvements  in  performance  were  most  likely  from  drag  reduction  with  useful  gains  from  bi-furcating  and 
minimizing  the  diameter  of  suspension  lines.  However,  the  greatest  gains  might  be  achieved  by  closing  the  leading 
edge  of  the  parachute  which  in  turn  necessitates  sweeping  the  wing  to  maintain  the  inflated  shape.  Glide  ratios  of  6:1 
were  predicted,  doubling  the  glide  ratio  of  existing  systems. 

Crimi  used  a  3 -DoF  model  (lateral  channel  only)  to  determine  the  relationship  of  aerodynamic  and  inertial 
parameters  to  lateral  stability  characteristics.28  He  represented  the  system  by  a  rigid  wing  with  a  suspended  payload 
in  a  steady  glide  and  analytically  prescribed  stability  boundary  for  spiral  divergence.  His  calculations  showed  that 
there  is  a  minimum  effective  anhedral  angle  required  for  absolute  stability.  Increasing  suspension  line  length  is 
stabilizing  for  spiral  divergence,  whereas  increasing  the  glide  slope  is  stabilizing  for  both  spiral  divergence  and 
oscillatory  response.  Iosilevskii  used  a  3 -DoF  model  (longitudinal  channel  only)  to  establish  center  of  gravity  and 
lift  coefficient  limits  for  a  gliding  personal  parachute.29 

Jann  discussed  two  different  mathematical  models  which  described  the  real  parafoil-payload  system.30  One  was 
a  3-DoF  model  and  another  was  a  4-DoF  model.  In  3-DoF  model,  only  the  motion  of  center  of  mass  is  considered. 
Also,  steady  state  longitudinal  motion  and  no  dynamic  effect  on  lateral  motion  were  considered.  The  4-DoF  model 
also  took  the  roll  angle  into  account. 

Machin  et  al.  developed  and  thoroughly  investigated  low-order  degree  of  freedom  models  for  longitudinal  and 
lateral  channels  (separately)  of  the  large-scale  ARS  (X-38)  for  the  space  station  crew  recovery  vehicle.22  They 
presented  a  method  of  extracting  system  longitudinal  and  lateral-directional  aero  data  to  feed  these  models.  These 
data/models  were  further  applied  to  a  NASA  Johnson  Space  Center  8-DoF  simulator  to  receive  a  reasonable 
dynamic  response  match  between  simulations  and  flight- test  data.31 

Prakash  considered  longitudinal  stability  of  the  parafoil  text  vehicle  ALEX  and  basically  used  one  degree  of 
freedom  model.32 

C.  Advanced  Modeling 

Being  a  very  light  structure  parafoil  is  strongly  influenced  by  the  apparent  mass  effects  of  the  surrounding 
air.8,33  The  majority  of  aforementioned  models  did  not  take  it  in  account  until  this  subject  was  discussed  by 
Lissaman  and  Brown,  who  modeled  the  parafoil  as  an  ellipsoid.34  Brown  provided  a  discussion  of  the  increasing 
effect  of  apparent  mass  on  parafoil  turning  performance  as  the  scale  gets  larger.35  Thomasson  derived  a  new 
formulation  of  the  EoM  of  a  rigid  body  in  a  moving  fluid  medium  which  might  have  velocity  gradients  or 
accelerations.36  He  separated  out  the  inertia,  apparent  mass,  and  relative  velocity  effect  from  derived  EoM.  Using 
these  EoM,  he  revealed  sources  of  error  in  existing  parafoil  model  using  classic  perfect  fluid  result.  Thus,  the 
formulation  provided  a  common  set  of  EoM  for  describing  the  motion  of  underwater  vehicles,  airships,  parafoils, 
and  airplanes. 
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Almost  all  studies  of  apparent  mass  effects  focused  on  an  ellipsoid,  a  two-dimensional  shape,  or  an 
axisymmetric  body.  An  ellipsoid  enjoys  the  virtue  of  having  three  planes  of  symmetry  and  a  geometric  centroid  that 
is  at  the  intersection  of  all  three  planes.  A  single  center  of  apparent  mass  can,  therefore,  be  found  that  is,  not 
surprisingly,  at  the  centroid.  Similarly,  when  there  is  axisymmetry,  the  center  of  apparent  mass  must  lie  along  the 
centerline,  and  it  turns  out  that  in  this  case  a  single  apparent  mass  center  can  also  be  found  that  is  at  the  volumetric 
centroid.  In  the  more  general  case,  it  may  not  be  possible  to  find  a  single  point  at  which  the  rotational  and 
translational  motions  are  decoupled. 

Barrows  examined  the  case  in  which  there  are  only  two  planes  of  symmetry,  such  that  the  parafoil  has  lateral 
symmetry  and  fore-and-aft  symmetry,  but  the  top  half  is  not  the  same  as  the  bottom  half.37,38  Therefore,  for  each  of 
the  three  axes  of  motion,  what  we  have  been  calling  the  center  must  lie  somewhere  on  the  intersection  of  the  two 
planes  of  symmetry,  but  they  may  not  coincide  allowing  for  more  than  one  center  of  apparent  mass.  Barrows 
computed  the  apparent  mass  of  a  parafoil  shape  for  motions  along  various  axes  using  potential  flow  analysis  (no 
account  has  been  taken  of  the  manner  in  which  the  presence  of  vortices  may  alter  the  kinetic  energy  of  the  flow). 
From  this,  the  6x6  apparent  mass  matrix  about  some  reference  point  was  derived  and  parametric  forms  for 
estimating  the  terms  of  this  matrix  were  suggested.  The  existence  of  multiple  mass  centers  results  in  off-diagonal 
terms  in  this  matrix  that  couple  the  translational  and  rotational  motions.  A  fully  populated  apparent  inertia  matrix 
has  36  terms,  of  which  only  21  are  distinct  because  it  must  be  symmetric.  If  the  body  has  both  lateral  and  fore-and- 
aft  symmetry,  there  are  only  eight  fundamental  parameters  of  apparent  inertia.  Barrows  has  shown  how  the  non¬ 
diagonal  6x6  apparent  inertia  matrix  about  a  certain  reference  point  can  be  used  to  compute  the  corresponding 
apparent  mass  matrix  at  any  other  reference  point.  He  also  presented  dynamic  EoM  including  nonlinear  terms. 

Barrows  noticed  that  whoever  stated  that  a  little  bit  of  knowledge  is  a  dangerous  thing  could  have  had  apparent 
mass  in  mind.  As  an  example  he  presents  comparison  of  parachute  roll  motions  about  a  point  near  the  confluence 
point  of  the  suspension  lines.  He  showed  that  the  effect  of  the  apparent  roll  inertia  may  be  greatly  reduced  by 
spanwise  camber.  Lissaman  and  Brown,  using  a  wing  model  with  no  spanwise  camber  (a  flat  wing),  concluded  that 
the  apparent  roll  inertia  may  be  as  much  as  five  times  the  solid  mass  roll  inertia.34  The  analysis  by  Barrows 
concluded  that  roll  inertia  of  a  cambered  wing  about  the  roll  center  may  be  as  small  as  0.1  of  the  corresponding 
value  for  a  flat  wing.  The  result  is  that  the  roll  acceleration  is  much  greater  than  expected  from  a  flat-wing  model. 
Steady-state  turn  rates  may  also  be  altered  when  compared  to  the  flat-wing  results,  but  not  as  dramatically. 

For  a  long  time  2-  DoF  or  3 -DoF  models  were  the  only  models  needed  to  analyze/optimize  the  integral  behavior 
and  design  of  parafoils.  Then,  the  need  for  a  higher  fidelity  models  emerged  to  study  the  behavior  parachute- 
submunition  systems  and  recently  to  support  autonomous  control  system  development.  Some  of  the  on  autonomous 
delivery  system  design  projects  are  the  Automatically  Guided  Parafoil,  by  European  Space  Agency  (1997),  Small 
Autonomous  Parafoil  Landing  Experiment  ALEX  by  Institute  of  Flight  Mechanics  of  the  German  Aerospace  Center 
(1999), 39,40  High  Altitude  Balloon  Experiments  in  Technology  HABET  by  Iowa  State  University,41  Automatic 
Parafoil  by  RIEGL  (2000),  Advanced  Precision  Aerial  Delivery  System  by  FXC  Corporation  (2001), 42  Precision 
Guided  Parachute  System  by  Atair  Aerospace  (2003)  and  others.  Obviously  for  more  sophisticated  models  the 
coefficients  provided  by  the  experimental  research  were  not  sufficient.  That’s  where  Prandtl  lifting  line  theory  and 
computational  fluid  dynamics  (CFD)  methods  helped  evaluate  missing  coefficients.43'48  For  instance,  Brown35  used 
scalar  aerodynamic  coefficients  calculated  using  discrete  vortex  Weissinger  method  (extended  lifting  line  theory)  by 
Kroo.49  (It  should  be  noted  that  while  the  coefficients  provided  by  wind-tunnel  experiments  are  know  within  the 
wide  range  of  angle  of  attack  and  side-slip  angle  the  coefficients  supplied  by  CFD  are  usually  just  static  scalars.) 

Detailed  dynamic  simulation  of  the  flight  mechanics  of  parachute  (or  parafoil*)  and  payload  systems  (two 
bodies  having  a  specific  joint  at  the  confluence  point)  appears  to  have  commenced  in  1968  with  PhD  work  of  Wolf, 
who  considered  the  stability  of  a  parachute  connected  to  a  load.50  Pillasch  and  Shen  concluded  that,  6-DoF  model, 
useful  for  rigid  body  dynamics,  is  inadequate  for  the  flight  dynamics  of  the  parachute- submunition  system,  therefore 
it  requires  more  degrees  of  freedom.51  They  reported  that  for  the  three-body  parachute  submunition  including  a 
friction  plate  and  elastic  suspension  lines  they  developed  two  versions  of  the  enhanced  model:  10-DoF  and  15-DoF 
(it  our  opinion  that  here  and  in  some  other  publications  degrees  of  freedom  for  multiple-body  systems  are  not  being 
counted  properly  but  we  leave  as  it  was  reported  in  the  reference  and  will  discuss  it  in  the  following  subsection). 
They  also  provided  Lagrangian  equations  and  some  details  on  the  first  one,  10-DoF  model.  In  1987  Reif  published 
another  9-DoF  model  and  simulation  results  which  compared  well  with  measured  oscillations  of  a  parachute- 
payload  system  that  was  mounted  on  a  van.52 


Once  a  6-DoF  model  of  parachute  or  parafoil  has  been  developed,  adding  two  or  three  DoF  (depending  on  a  specific  joint)  for 
payload  is  a  relatively  minor  addition,  which  is  quite  similar  for  both  parachutes  or  parafoils.  More  enhanced  model  may  include 
the  model  of  a  joint  (say  elastic  riser)  itself  allowing  for  more  DoF. 
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The  parachute-riser-payload  system  by  Fallon  was  modeled  as  three  bodies  each  with  6-DoF.53  Since  the 
parachute  and  payload  were  connected  by  the  riser,  and  given  the  constraints  of  the  riser,  they  reduced  parachute- 
riser-payload  system  to  a  15-DoF  problem.  Later,  Doherr  and  Schilling  reported  on  the  development  of  a  9-DoF 
dynamic  model  of  rotating  parachutes  which  was  used  to  decelerate  and  spin  the  submunitions  of  defense 
weapons.54,55  This  model  described  the  parachute’s  behavior  using  Euler  differential  equations  and  payload’s 
attitude  using  quaternions  (to  avoid  singularities).  It  therefore  included  19  (versus  usual  18)  independent  variables. 
The  differential  equations  were  later  converted  to  state-space  representation  by  considering  the  apparent  mass  terms 
to  be  external  forces  and  moments  caused  by  accelerations  and  velocities  at  an  earlier  time  (i t-z )  (the  time  delay  r  was 
derived  from  the  earlier  experimental  observations  during  water  tunnel  tests).  The  stochastic  disturbance  was  added 
to  aerodynamic  forces  to  address  the  issue  of  nonrigid  nature  of  the  parachute  system.  By  comparing  results  from  6- 
DoF  and  9-DoF  models,  they  conclude  a  9-DoF  model  predicts  stability  characteristics  more  adequately 
(realistically).  Furthermore,  their  work  established  the  very  high  sensitivity  of  the  motion  of  a  parachute-payload 
system  to  atmospheric  disturbances. 

Cockrell  and  Haidar  considered  the  influence  of  canopy-payload  coupling  (accounting  for  added  masses)  on 
dynamic  stability  in  pitch.56  They  used  planar  equations  for  two  bodies  (they  refer  to  the  earlier  work  by  Frucht57). 
The  paper  by  Moulin  pointed  out  the  importance  of  the  modeling  the  link  between  the  parachute  and  the  load  by 
showing  its  contribution  to  the  dynamic  behavior  of  this  complex  system.58  An  8-DoF  model  was  presented  with  its 
sensitivity  to  the  driving  parameters  based  on  physical  considerations.  Vishnyak  considered  12-DoF  model  of  non- 
rigid  parachute-payload  system  with  elastic  connection  (they  used  quaternions  for  kinematics).59 

Hailiang  and  Zizeng  used  a  9-DoF  model  to  study  motion  of  a  parafoil  and  payload  system.60  In  contrast  to 
Doherr  and  Schilling,  they  reported  only  small  differences  in  the  motion  and  stability  between  6-DoF  and  9-DoF 
dynamic  models.  In  studying  stability  characteristics  as  a  function  of  the  pitch  inertia  of  the  pay-load,  Hailiang  and 
Zizeng  found  the  decay  ratio  and  period  increase  as  pitch  inertia  is  increased. 

The  latest  model  of  Wolf  allowed  5  DoF  for  each  body  of  the  parachute-payload  system,  with  roll  about 
axisymmetry  being  neglected.61  The  major  assumptions  are  rigid,  axisymmetric  parachute  and  payload  bodies,  flat 
earth  with  no  wind.  The  model  described  motion  of  two  dynamically  coupled  rigid  bodies  upon  which  simplified 
aerodynamic  forces  acted.  Using  this  10-DoF  representation,  Wolf  established  that  stability  is  reduced  as  riser  length 
is  increased  or  parachute  weight  is  increased  and  that  stability  is  improved  by  increasing  parachute  axial  and  normal 
aerodynamic  force.  More  recent  efforts  by  Zhu  et  al.  as  well  as  Gupta  et  al.  have  incorporated  parafoil  structural 
dynamics  into  the  dynamic  model  of  a  parachute  and  pay-load  system.47,62 

In  support  of  future  control  system  development  Entchev  et  al.  characterized  the  dynamics  of  a  small  parafoil 
(canopy  and  payload)  as  a  6-DoF  nonlinear  numerical  simulation  based  on  a  rigid-body  paradigm  (aircraft  EoM)  and 
aerodynamic  model  of  the  system  (23  static  coefficients).63  This  simulation  was  based  on  construction  of  a  scale 
rigid  model  of  the  canopy  and  wind  tunnel  testing  of  the  same,  CFD  analysis,  vehicle  geometry  measurements  and 
mass  properties  estimation.  The  validity  of  the  simulation  was  assessed  using  further  CFD  analysis  and  comparison 
with  current  parafoil  theory.  Sensitivity  to  various  parameters  was  established  and  the  accuracy  of  the  rigid  body 
approximation  was  assessed  as  quite  satisfactory  for  further  design  of  guidance,  navigation  and  control  (GNC) 
system  for  guided  Airborne  Delivery  System  (GPADS  for  the  US  Army  Soldier  System  Command).  Mortaloni  et  al. 
reported  the  development  of  a  6-DoF  Simulink  model  of  the  Pegasus-500  parafoil-payload  system  by  FXC 
Corporation. 10  Both  diagonal  and  non-diagonal  terms  of  the  apparent  mass  matrix  were  taken  into  account.  Model 
was  tuned  to  perfectly  match  results  of  a  real  drop  of  the  system.  The  complete  model  was  later  used  for  the  design 
and  simulation  of  the  autonomous  GNC  system. 

Siegers  and  Costello  presented  a  model  in  which  combined  system  of  parafoil  and  payload  was  modeled  with  9 
DoF.64  The  model  included  three  inertial  position  components  of  the  joint  (incorporating  rotational  stiffness  and 
damping  coefficients  of  joint)  as  well  as  the  three  Euler  orientation  angles  for  parafoil  and  payload.  The  model  used 
apparent  terms  of  Lissaman  and  Brown.34  The  parafoiFs  shape  was  modeled  as  a  collection  of  joint  panels  (Fig.  1)  to 
address  the  issue  of  realistic  initial  disturbance  while  X-38  parafoil  section-by-section  opening.  It  was  shown  that  for 
the  small  brake  deflections  parafoil-payload  systems  exhibit  two  modes  of  directional  control:  roll  and  skid  steering. 
Although  the  model  only  employed  basic  aerodynamic  coefficients,  neglecting  the  majority  of  them,  it  could 
however  simulate  reasonable  control  response  from  both  left  and  right  deflection  and  parafoil  tilt.  Using  a  dynamic 
modeling  approach  similar  to  Doherr  and  Schilling  and  Hailiang  and  Zizeng,  stability  and  control  characteristics  of 
this  system  were  examined.  Particular  attention  was  paid  to  steady-state  control  response  as  a  function  of 
fundamental  design  parameters  such  as  parafoil  canopy  geometry,  angle  of  incidence,  and  varying  control 
deflection. 

Mooij  et  al.  described  detailed  mathematical  modeling  of  9-DoF  parafoil  and  payload  system.65  The  parafoil 
system  was  modeled  as  a  rigid  parafoil  and  rigid  store  that  were  connected  to  each  other  by  rigid  suspension  lines 
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and  a  riser.  Included  in  the  riser  was  a  hinge  that  allowed  free  three  axes  rotation  of  the  parafoil  with  respect  to  the 
store  and  it  transfers  only  forces,  but  no  moment.  The  hinge  was  modeled  as  a  damped  spring.  The  earth  was 
assumed  to  be  spherical  and  rotating  to  account  for  Coriolis  and  transportation  acceleration. 


Figure  1.  Parafoil  canopy  geometry  used  by  Siegers  and  Costello.64 


The  latest  work  by  Muller  et  al  described  the  highly-realistic  8 -DoF  model  of  the  parafoil  system.66  It  presented 
nonlinear  differential  equations  of  the  system  and  included  detailed  simulations  of  dynamic  reactions  to  symmetric 
and  asymmetric  control  inputs,  to  atmospheric  disturbances  and  to  a  twisted  payload. 


III.  Identification  and  Scalability 

As  seen  from  the  previous  section,  a  variety  of  models  of  different  parafoil-payload  systems  has  been  developed 
over  past  40  years.  The  majority  of  them  deals  with  simplified  models  addressing  some  specific  issues.  However, 
several  higher- fidelity  models  have  also  been  presented.  By  higher-fidelity  we  mean  6-DoF  models  considering 
parafoil-payload  system  as  a  rigid  body,  and  higher  DoF  models  considering  a  flexible  joint  between  the  parafoil 
and  payload.  Some  of  these  joints  allow  pitching,  rolling  and  yawing  of  the  payload  around  the  confluence  point 
(joint)  therefore  adding  3  additional  DoF.  So  in  this  case  9-DoF  models  are  considered  (alternatively  one  may  say  it 
is  a  6-DoF  model  for  the  payload  plus  3  DoF  for  parafoil  with  respect  to  the  payload).  If  the  joint  restricts  some  of 
these  motions  (usually  rolling)  then  it  is  lower  than  9-DoF  models  (8-DoF  for  non-rolling  joint).  Introducing  a  joint 
or  a  suspension  system  (for  instance  having  elastic  suspension  lines)  as  a  separate  subsystem  obviously  adds  more 
DoF. 

Some  of  the  aerodynamic  coefficients  in  these  high-fidelity  models  are  taken/extrapolated  from  aforementioned 
experimental  wind-tunnel  data  (meaning  they  do  not  belong  to  the  concrete  system  one  tries  to  model)  and/or  CFD 
simulations  (which  are  not  complete  and  may  be  way  off),  some  are  simply  neglected.  Some  papers  provide  results 
of  flight  tests  and  suggest  how  to  extract  and  validate  certain  characteristics.  The  issue  of  apparent  masses  has  also 
been  addressed  and  empirical  formulas  for  both  diagonal  and  non-diagonal  terms  have  been  derived  and  used  by 
many  researchers  in  their  models. 

The  models  developed  for  different  parafoil-payload  systems  behave  well  and  produce  reasonable  transition 
behavior  and  integral  performance,  attesting  quite  a  mature  stage  of  such  systems  modeling.  Therefore,  one  of  the 
important  issues  to  address  now  is  not  the  development  of  the  model  itself  (besides,  EoM  are  well-known  and  initial 
guesses  on  the  majority  of  the  coefficients  can  be  also  obtained),  but  the  technique  capable  of  reliable  and  effective 
model  identification  with  consecutive  adjustment  of  initial-guess  coefficients  to  better  approximate  real  drop  data. 
Recently,  the  scaling  a  parafoil-payload  model  (its  coefficients  and  parameters)  also  became  an  important  issue.  (To 
our  knowledge  only  Goodrick  presented  data  and  analysis  of  scale  effects  on  ram-air  parachutes  based  on  his  6-DoF 
model  simulations  for  parafoil  wings  areas  from  32  to  6400 ft2  and  for  wing  loadings  from  1  to  4 lb/ft2.61) 

As  mentioned  in  the  previous  section  the  variety  of  models  higher  than  the  point-mass  models  exists  for  both 
parafoil-payload  systems: 

•  9-DoF  model  by  Hailiang  and  Zizeng;60 

•  8-DoF  model  by  Siegers  and  Costello;64 

•  9-DoF  model  by  Mooij  et  al.;65 

•  8-DoF  model  by  Muller  et  al.  ;66 

and  parachute-payload  systems  (the  order  of  the  model  corresponds  to  that  of  claimed  by  the  authors): 

•  10-DoF  model  by  Pillasch  and  Shen;51 

•  9-DoF  model  by  Reif;52 

•  9-DoF  model  by  Doherr  and  Schilling;54,55 

•  8-DoF  model  by  Moulin;58 

•  12-DoF  model  by  Vishnyak;59 

•  15 -DoF  model  by  Fallon.53 
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Moreover,  if  the  6-DoF  model  of  some  system  has  been  developed  already,  than  adding  additional  degrees  of 
freedom  for  the  payload  is  a  natural  next  step  (the  problem  is  not  in  adding  two  or  three  more  degrees  rather  than  in 
identification  of  such  model  since  usually  the  only  data  recorded  during  the  real  drop  is  the  data  related  to  the 
payload).  The  EoM  written  for  the  center  of  gravity  of  the  entire  system  (Fig.2a),  should  be  rewritten  for  the  parafoil 
and  payload  separately  (Fig.2b)  accommodating  additional  equations  describing  parafoil-payload  coupling  reactions 
for  the  concrete  bundle  (at  the  confluence  point).  So  the  previously  developed  and  well-tested  6-DoF  model  of  the 
Pegasus  parafoil-payload  system10  was  used  as  a  basis  for  the  scalable  model  development. 


Figure  2.  Rigid  one-body  system  (a)  versus  flexible  two-body  system  (b).18 


In  general  the  8-DoF  (9-DoF)  model  of  the  parafoil-payload  system  can  be  brought  to  the  following  sets 
dynamic  and  kinematic  equations 


V*r  =  A”1 


r  f  ’ 

-LAV* 

K_ 

(1) 


At=Re(A>c 
Ap  =R£(Ap)mp 

P =  *RV , 


(2) 

(3) 

(4) 


where  vector  V*  =  J  consists  of  a  velocity  vector  of  the  confluence  point  V  =  [u,v,wj ,  angular 

velocity  vectors  of  parafoil  (canopy)  (ac  =[pc,qc,rcJ  and  payload  (ap  =  [o,qp,rp J  ,  vectors  Ac  and  Ac 


(A  =  \q>,9,y/f )  constitute  Euler  angle  triads  for  the  canopy  and  payload  respectively,  P;  presents  an  inertial 


position  of  the  system  (confluence  point).  On  the  right-hand  side  of  Eq.(l),  block  3-by-3  matrix  A  (where  each 
element  is  a  3-by-3  matrix  itself)  represents  a  mass-inertia  matrix  of  the  system,  vectors  F,  Mc  and  constitute 
aerodynamic  force  and  moment  acting  on  the  system  and  its  two  components  (canopy  and  payload  respectively),  and 
matrix  £  is  a  block  3-by-3  matrix  where  each  element  is  a  3-by-3  matrix  too.  In  Eqs.(2),(3),  notation  R^(A)  stands 
for  the  matrix  operator  acting  on  a  vector  A 


A  = 


1 

0 

0 


sin  cptgO 
cos  (p 

sin  cp!  cos  6 


cos^tg^ 
-sin  cp 
cos  cp!  cos  6 


CO  =  R£(A)« 


(5) 


and  matrix  in  Eq.(4)  is  the  rotation  matrix 
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iR  = 


cos  yj  cos  6 

cos  yj  sin  6  sin  cp-  sin  yr  cos  cp 
cos  yf  sin  6  cos  cp  +  sin  y/  sin  cp 


sin  yj  cos  6 

sin  yj  sin  6  sin  cp  +  cos  yj  cos  cp 
sin  yj  sin  6  cos  cp-  cos  yj  sin  cp 


-sin(9 
cos  6  sin  cp 
cos  6  cos  cp 


(6) 


Equation  (1)  (nine  scalar  equations)  describes  systems  dynamics,  while  the  rest  of  them  (Eqs.(2)-(4))  -  system’s 
kinematics.  If  the  opening  stage  is  to  be  also  considered  than  in  order  to  avoid  singularities  it  is  worth  substituting 
Eq.(3)  with 


Q, 


0 

CO 


(7) 


(notation  S(£)  = 
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£ 
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A 
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o 


for  the  skew  matrices  is  used),  applicable  for  quaternions 
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sin— ^ cos  — cos - 
2  2  2 


-cos  —  sin  —  sin- 
2  2  2 


=  Rp(AJ 


(8) 


(notation  R0(A)  stands  for  the  matrix  operator  acting  on  a  vector  A ).  The  rotation  matrix  in  Eq.(4)  in  this 
case  would  be  compound  as 


1  —  2  {ql+ql) 


;r  = 


2(<M,  -?o^3) 
2{q,qx+qaq2) 


2(qtq2+qltqi) 

l-2(q23+qf) 

2(?3?2-?0?l) 


2(?,?3-?0?2) 
2(q2q3+q0ql)  , 

1  -  2(q l  +ql) 


(9) 


and  the  relationship  between  quaternions  and  Euler  angles  established  by 


A 


p 


a,,-,  22miMiL 

1-2  (ql+ql) 


-sin_1(2  (qiq3-q0q2)) 


tan-l  2(glg2+gQ <h) 
l-2(ql+q]) 


(10) 


The  dependences  of  aerodynamic  and  control  coefficients  versus  angle  of  attack  were  developed  earlier  based 
on  the  available  data13  tuned  to  match  the  flight  test  data.10  Parafoil-payload  coupling  reactions  were  modeled  in 
pretty  much  the  same  manner  as  say  in  Refs. 54, 55, 66.  (As  expected  at  this  time,  not  many  changes  were  observed 
when  the  outputs  of  the  6-DoF  and  8 -DoF  models  were  compared.) 

Further  attempts  were  devoted  to  the  development  of  a  scalable  model  assuming  that  the  data  from  real  drops  of 
several  different  parafoil-payload  systems  becomes  available.  Although  up  until  now  only  two  incomplete  sets  of 
data  for  each  of  two  systems  shown  on  Fig. 3  is  available  (both  systems  only  differ  by  the  way  payload  is  attached) 
the  following  addresses  some  steps  made  to  convert  the  model  of  a  specific  payload  delivery  system  to  that  of 
general  system  requiring  only  limited  amount  of  geometry/mass  inputs. 


Figure  3.  Generic  parafoil-payload  system  with  different  ways  of  payload  attachment. 
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We  assume  first  that  nothing  but  just  a  single  geometric  parameter  is  known  a  priory.  The  goal  then  would  be  to 
determine  the  rest  of  them  from  observation  of  the  descending  system.  To  this  end,  Fig.4  shows  an  example  of 
several  front- view  images  of  the  system  with  superimposed  theoretical  skeleton. 


Figure  4.  Front- view  geometry. 


For  this  skeleton  the  following  theoretical  dependences  hold 


R  = (2o0_ 


b  = 


2{R- 


a  jtga  . 


b2  =SR2(l-cosa)-4a2 


(11) 


(here  a  bar  denotes  relative  quantities  related  to  the  span  of  uninflated  parafoil  bu).  These  dependences  of  front- view 
system  parameters  versus  flare  angle  are  shown  on  Fig. 5.  The  values  of  the  very  same  parameters  as  estimated  from 
analyzing  video  data  are  also  presented  on  Fig. 5  demonstrating  the  perfect  match.  That  means  that  by  observing  the 
flare  angle  with  known  span  bu  the  rest  parameters  (radius  R ,  arch  a ,  span  b ,  arch-to-span  ratio  and  front  suspension 
line  length)  can  be  accurately  estimated  using  Eqs.(l  1).  The  accuracy  of  the  estimates  can  be  found  from  Eq.(l  1)  as 

SR=_Sa  (12) 

R  a 


Figure  5.  Dependences  of  relative  front-view  geometry  parameters  on  the  flare  angle. 


Proceeding  with  analyzing  the  side-view  geometry  of  the  system  with  now  known  radius  R ,  it  is  possible  to 
reconstruct  the  magnitude  of  a  chord  c,  aspect  ratio  b/c  and  the  geometry  of  the  suspension  system  (Fig.6).  As 
known  the  aspect  ratio  defines  the  parafoil  aerodynamics,  so  that  the  nominal  dependences  of  aerodynamic  and 
control  coefficients  versus  angle  of  attack  (developed  originally  for  the  certain  system  with  a  certain  aspect  ratio) 
may  now  be  adjusted  by  some  empirical  coefficients  to  be  found  during  the  following  system  identification  stage. 

Knowing  the  front-  and  side-view  geometry  also  allows  to  proceed  with  the  defining  a  location  of  the  center  of 
gravity  (CG)  of  the  airfoil  and  inflated  parafoil  (assuming  the  circular  model,  center  of  gravity  of  the  inflated 
i  sin  cc  i 

parafoil  lies  R  1 - beneath  the  CG  of  the  airfoil  along  the  line  tilted  clockwise  by  the  rigging  angle)  as  well 

\  a  ) 

with  the  estimating  parafoil  and  entrapped  air  masses  (Fig. 7). 

The  following  step  would  be  adding  payload  data  and  defining  the  type  of  a  joint  to  adjust  coupling  reactions 
between  parafoil  and  payload. 
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Figure  6.  Side-view  geometry  (a)  and  geometry  of  the  suspension  system  (b). 


Figure  7.  Side-view  geometry  of  the  airfoil. 

The  set  of  variable  gains  was  introduced  to  address  the  model  identification/scalability  issue.  This  set  includes  a 
group  of  aerodynamic  force  coefficients  kF  ,  a  group  of  moment  coefficients  kM ,  a  group  of  mass-geometry 
coefficients  kG ,  and  a  group  of  apparent  mass  coefficients  ka .  The  total  number  of  these  coefficients  exceeds 

several  dozens.  Nominally,  for  the  original  model  these  coefficients  or  gains  are  equal  to  unity.  The  goal  of  the 
following  identification  is  to  adjust  these  coefficients  so  that  they  provide  the  best  match  between  simulations  and 
real  drop  data  for  several  different  parafoil-payload  systems  simultaneously.  In  such  a  case  the  defined  gains  would 
constitute  a  generic  parafoil-payload  system  model  allowing  to  use  it  as  an  initial  guess  when  designing  a  brand  new 
system  (when  only  a  few  basic  parameters  are  known).  Of  course  later  the  model  may  be  adjusted  further  to  match 
the  real  drop  data  of  this  specific  system  rather  than  some  generalized  system. 

Figure  8  presents  the  hierarchy  and  the  task  flow  of  the  developed  software. 


SingleRun.m 


Choose  Data  Set 
Choose  Scale  Factor 


Geometry.m 


Define  Geometry 


OptimalGains.m  ) 

Interactive  Interface 
Investigate  Configuration 


Spectrums.m _ | 

Define  Multiple  Criteria 


LoadVariables.m 


LoadReferenceData.m 


LoadAtmosphere.m 


LoadlnitialConditions.m 


SystemParameters.m 


NominalParameters.m 


GetGains.m 


CenterofGravity.m 


Masslnertia.m 


ADSModel.mdl 


Errors.m 


Graphs.m 


MinSearch.m 


DeflneGlobalVariables.m 


Single  Criterion  Routines 
Multiple  Criteria  Identification 


Figure  8.  Software  architecture. 
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The  core  element  in  this  hierarchy  is  the  Simulink  model  of  the  system  itself  ( ADSMODEL.mdl  on  Fig.8).  But 
before  this  system  can  be  run,  all  input  parameters  should  be  supplied  (by  default  the  nominal  Pegasus  system  is 
considered).  If  the  real  drop  data  is  available  than  identification  routine  may  be  employed.  In  more  general  case 
when  almost  nothing  is  known  yet  for  sure,  multiparameter  multicriteria  software  is  used.11  If  there  is  a  need  to 
adjust  just  a  few  parameters  rather  than  dozens  of  them,  than  a  simple  single-criterion  method  can  be  employed. 
Figure  9  shows  that  during  identification  (or  after  it)  any  variable  gain  may  be  interactively  changed  directly  that 
gives  the  developed  tool  sort  of  flexibility  from  one  the  hand  and  the  way  to  explore  the  influence  of  a  certain 
parameter  from  the  other  hand. 


IV.  Conclusion 

The  paper  describes  the  previous  efforts  in  modeling  different  parafoil-based  payload  delivery  systems  and  in 
obtaining  aerodynamic  and  control  coefficients  for  the  parafoil  itself.  Based  on  the  available  wind-tunnel  data  for 
different  parafoils  and  recent  successes  in  the  development  and  validation  of  the  high-fidelity  models  of  several 
similar  parafoil-payload  systems,  the  efforts  have  been  made  to  convert  the  6-DoF  model  developed  earlier  for  a 
specific  parafoil-payload  system  to  a  higher-degree-of-freedom  model  capable  to  be  tuned  to  another  system, 
meaning  capable  to  be  scaled.  Several  dozens  of  variable  gains  were  introduced  to  allow  this  capability.  While  the 
paper  presents  a  current  status  of  development  of  correspondent  software  to  address  the  issue  of  model  identification 
and  scalability,  the  author  hopes  to  use  this  software  to  proceed  with  a  scalable  high-fidelity  model  of  the  generic 
parafoil-payload  system  in  the  nearest  future  when  more  experimental  data  becomes  available. 
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